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Abstract 

A new method for analyzing complex superhump light curves for the 2001 outburst of WZ Sagittae 
is proposed. The complexity arises because intrinsically time-varying and non-axisymmetric distributions 
of superhump light sources are coupled with the aspect effects around the binary orbital phase because 
of its high orbital inclination. The new method can disentangle these complexities by separating the 
non-axisymmctric spatial distribution in the disk from the time variation with the superhump period. It 
may be called a helical tomography of an accretion disk because it can reconstruct a series of disk images 
(i.e., disk's azimuthal structures) at different superhump phases. The power spectral data of superhump 
hght curves of the 2001 outburst of WZ Sge by Patterson et al.(2002,PASP,114,721) are now interpreted 
under a new light based on the concept of helical tomography, and the azimuthal wave numbers of various 
frequency modes are identified. In particular, a frequency component, nujQ — 51, where cjq and Q. are the 
orbital frequency and a low frequency of the apsidal precession of the eccentric disk, is understood as an 
(n — l)-armed traveling wave in the disk. A vigorous excitation of a wave component of cos(20 — Swo^) in 
the first week of the superhump era of WZ Sge, where 8 is the azimuthal angle, supports Lubow's (1991) 
theory of non-linear wave coupling of the eccentric Lindblad resonance for the superhump phenomenon. 
This method can in principle be applied to other SU UMa stars with high orbital inclination if light curves 
are fully covered over the beat cycle. 

Key words: accretion, accretion disks — stars: binaries: close — stars: dwarf novae — stars: individual 
(WZ Sge) — stars: nova, cataclysmic variables 



1. Introduction 

WZ Sagittae, a prototype of one subclass of dwarf no- 
vae, called WZ Sge, underwent a full-scale outburst in 
2001 July and extensive observations, particularly in op- 
tical light, covering almost all nights without any signif- 
icant gaps, were made for more than one hundred days 
from the start of eruption on July 23 to the long fading 
to quiescence with a worldwide collaboration of amateur 
and professional astronomers. Its beautiful light curves 
can be found in VSNET (htpp://www.kusastro.kyoto- 
u. ac.jp/vsnet.DNe/wzsge01. html) and AAVSO, and some 
of them have already been published by Ishioka et al. 
(2002) and by Patterson et al (2002). The 2001 super- 
outburst of WZ Sge exhibited many interesting features, 
which will surely be discussed in the future by various 
workers. This paper addresses only one of those features, 
i.e., an analysis of complex light curves of common super- 
humps observed during the main outburst. 

WZ Sge exhibited in its photometric light curves peri- 
odic humps repeating with a binary orbital period called 
either "outburst orbital humps" by Patterson et al (2002) 
or "early superhumps" by Ishioka et al (2002) or "early 
humps" by Osaki and Meyer (2002) for the first twelve 
days in its 2001 outburst. The origin of this phenomenon 
was discussed by Osaki and Meyer (2002) and by Kato 
(2002), and will not be discussed in this paper. 



After the first twelve days, the common superhumps 
with a period of 1% longer than the binary orbit emerged 
in accordance with the SU UMa-type nature of WZ Sge 
stars. The superhumps then subsisted all the way to 
the end of the main outburst and during the so-called 
echo outbursts and in the long fading stage to quiescence 
(Patterson et al 2002), although they may be of "late su- 
perhump nature" in its later stage. The overall picture of 
the 2001 outburst of WZ Sge based on the disk instability 
model (i.e., the thermal-tidal instability model) was pre- 
sented by Osaki and Meyer (2003) and may not need to 
be repeated here. 

Once superhumps emerge, light curves of WZ Sge be- 
come extremely complex, as seen, e.g., in figure 5 of 
Patterson et al. (2002). Although they look extremely 
complex, one remarkable fact may be recognized, that is, 
the basic patterns of light curves were more or less re- 
peated with the beat cycle of about six days between the 
superhump and the orbital periods. The reason why the 
superhump light curves in WZ Sge are so very complex is 
that intrinsically time- varying and non-axisymmetric dis- 
tributions of superhump light sources are coupled with the 
aspect effects around the binary orbital phases because of 
its high orbital inclination (Osaki, Meyer 2003). 

In this paper we propose a new method of light curve 
analysis to disentangle these complexities by separating 
the non-axisymmetric spatial distribution of the disk from 



2 



Y. Osaki 



a time variation occurring with the superhump period. In 
section 2, we explain the basic concept of helical tomog- 
raphy and how to reconstruct an "orbital light curve" at 
a given superhump phase. In section 3, we present the 
inversion method of an "orbital light curve" at a given su- 
perhump phase into an azimuthal brightness distribution 
of the disk. In section 4 we interpret the power spec- 
tra of light curves in the common superhump era of the 
2001 outburst of WZ Sge observed by Patterson et al. 
(2002) under a new light based on the concept of helical 
tomography, and we identify the azimuthal wave numbers 
of various frequency modes. In section 5 the modes ob- 
served in WZ Sge are interpreted in the context of Lubow's 
(1991) theory of non-linear wave interaction for the super- 
hump phenomenon. Possible applications of the method 
to other cataclysmic variable stars are dealt with in sec- 
tion 6. Section 7 is a summary and discussions. 

2. Basic Concept of Helical Tomography 

It is now well established that the superhump phe- 
nomenon is produced by a tidally driven eccentric insta- 
bility within accretion disks in cataclysmic variable stars 
(Whitehurst 1988; Hirose, Osaki 1990; Lubow 1991a,b). 
This instability occurs when an accretion disk expands to 
reach the 3:1 resonance region that excites the precessing 
eccentric structure in the accretion disk. 

We first note that the superhump phenomenon is pro- 
duced by the tidal stressing of a precessing eccentric disk 
by an orbiting secondary star, and that the intrinsic varia- 
tion within the accretion disk is strictly periodic with the 
superhump period. Wc may see how light distributions in 
the disk of the superhump binary system vary with the 
superhump phase in hydrodynamic simulations [c. g., in 
figure 7b of Hirose and Osaki (1990), figure 8 of Murray 
(1998)]. In fact, observed superhump light variations in 
many SU UMa stars (those of low orbital inclination) 
are singly periodic. However, a complication arises when 
the superhump phenomenon occurs in cataclysmic vari- 
ables with high orbital inclination because aspect effects 
cause other light variations in a non-axisymmctric disk, as 
shown by Simpson, Wood, and Burke (1998). These com- 
plications can be solved by the following method, which 
we call "helical tomography" . 

The method, itself, is very simple in principle. We re- 
produce at the first step an "orbital light curve" for a cer- 
tain superhump phase, and then invert the "orbital light 
curve" thus obtained into brightness distribution in the 
disk for the given superhump phase. In this section wc 
explain the method for reconstructing an "orbital light 
curve" for a given superhump phase. The inversion to an 
azimuthal brightness distribution in the disk is described 
in the next section. 

In order to reconstruct an "orbital light curve" for a 
given superhump phase, wc use all of the superhump light 
curves for one beat cycle of the orbital and superhump 
periods (i.e., about 6 days in the particular case of WZ 
Sge). Figure 1 illustrates the basic process of this method. 
Let us consider an ideal case in which superhump light 
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Fig. 1. Schematic illustrations showing how the helical to- 
mography works. The cylinder in the left stands for light 
curve observations over one superhump beat cycle. The cir- 
cle of the cylinder corresponds to one orbital cycle, while the 
beat phase proceeds from the top to the bottom of the cylin- 
der. Let the number of orbital light curves in one beat cycle 
be N. The cylinder is then a sum of N thin cylinders. Each 
thin cylinder is thought to show one orbital light curve at a 
certain beat phase, although light curves are not shown here. 
In our helical tomography, we cut a small piece from each 
orbital light curve at a given superhump phase, 'PsK^ pro- 
ceeding from the top to the bottom as illustrated by a helix 
in the left figure. We then assemble them into a fictitious 
"orbital light curve" at a given superhump phase as shown 
in the right figure. In this example we get ten "orbital light 
curves" for ten superhump phases. Each "orbital light curve" 
is then inverted to azimuthal light distribution for a given su- 
perhump phase by a method described in section 3. We call 
this method "helical tomography" because it can reconstruct 
a series of disk images (i.e., azimuthal light distribution) at 
different superhump phases, and because we cut out pieces 
of light curves "helically" to construct new light curves for 
certain superhump phases. 

curves are completely covered over the superhump beat 
cycle, i.e., in the case of WZ Sge, with an orbital period 
of 0.05669 d and a superhump beat period of about 6 d, 
which results in about 100 orbital light curves. We then 
align these one hundred light curves with increasing su- 
perhump beat phase, just like figure 5 of Patterson et al. 
(2002). We then cut them into small pieces in which each 
piece corresponds to a certain superhump phase. For in- 
stance, let us consider to reproduce ten light curves for 
ten superhump phases, ipsn- We then cut out ten pieces 
from each light curve, assemble them for a certain super- 
hump phase with increasing beat phase, and paste them 
as a function of the orbital phase to construct new light 
curves (in this case ten light curves). 

Each newly formed light curve is a fictitious "orbital 
light curve" at a given superhump phase, and can be in- 
verted into an azimuthal brightness distribution of the 
disk by a method described in the next section. This then 
gives a disk image at a given superhump phase seen from a 
particular inclination of the binary. The advantage of this 
method is that we can disentangle the geometrical effects 
of the azimuthal light distribution for a given superhump 
phase from the intrinsic time variation of the superhump. 
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itself. In other words, we translate the time variations of 
original light-curves into azimuthal light distributions at 
certain supcrhump phases. We call this method "helical 
tomography" of an accretion disk, because it can recon- 
struct a series of disk images (i.e., azimuthal light dis- 
tributions) at different superhump phases. The "helical 
tomography" is a medical terminology by which a three- 
dimensional image of a human body is reconstructed. We 
use this medical terminology because the method can be 
used to reconstruct a series of disk images in different 
superhump phases. The adjective "helical" is used in an- 
other sense in that we cut out "helically" pieces of light 
curves to construct new light curves for certain supcrhump 
phases (see, figure 1). 

3. Inversion of an "Orbital Light Curve" at 
a Given Superhump Phase into Azimuthal 
Brightness Distribution of the Disk 

3.1. A Simple Model for Superhump Light Distribution 

In order to extract information about the supcrhump 
light distribution in the accretion disk from an "orbital 
light curve" at a given superhump phase obtained by the 
method described in the preceding section, we need to 
make some simplifying assumption concerning the super- 
hump light source. Radiation from the accretion disk of 
SU UMa-type dwarf novae during the common superhump 
era consists of two different types of radiation. One is ra- 
diation from the surface of the whole accretion disk pro- 
duced by the ordinary viscous dissipation; the other is 
that from tidal dissipation near the disk edge. When we 
discuss superhump light curves, it is the common prac- 
tice that we are interested only in time-varying compo- 
nent with a time-scale less than the orbital period and 
we remove from the analysis secular time variation with 
a time-scale longer than the orbital period. In such cases 
we assume that the time- varying part of the radiation is 
solely produced by the tidal dissipation of the eccentric 
precessing disk, and that radiation from the main part of 
the accretion disk due to the ordinary viscous dissipation 
contributes only to the constant background radiation in 
superhump light curves. 

To extract information about the spatial distribution of 
superhump light at a given superhump phase from "orbital 
light curves" , we need to make further assumptions. Since 
an orbital light curve includes only "one-dimensional" in- 
formation, we can extract only one-dimensional informa- 
tion of the disk. Admittedly, the terminology "helical to- 
mography" may be an over-statement in this sense be- 
cause we are not able to reconstruct two-dimensional im- 
ages of the disk. However, we keep using this terminology 
because we can reconstruct a series of the disk images with 
different superhump phases. We adopt a simplest model 
for radiation from the tidal dissipation in a given super- 
hump phase. As discussed by Osaki and Meyer (2003), 
the tidal dissipation at a given superhump phase is non- 
axisymmetric in its pattern, and occurs near the disk rim 
in a region of about 1/10 width in the disk radius, i.e., al- 
most the same size as the estimated half thickness of the 
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disk at its rim. Thus, about one half of the superhump 
light may be expected to be radiated from the two hor- 
izontal disk surfaces and one half from the vertical disk 
rim. For binary systems with a high orbital inclination, 
radiation from the tidal dissipation comes mostly from the 
vertical disk rim; its relative importance of the disk rim to 
the horizontal disk surface is proportional to tani, where 
i is the orbital inclination. Furthermore, we assume for 
simplicity that the surface of the accretion disk is flat, so 
that radiation from the horizontal disk surface for a given 
superhump phase is independent of the orbital aspect. In 
other words, we may consider only superhump light from 
the vertical disk rim and we may regard radiation from 
the disk surface as the constant background light with the 
above-mentioned assumption when we discuss its aspect 
dependence around the binary orbit. 

3.2. Forward Problem: Transformation of Spatial 
Distribution of the Superhump Light at a Given 
Superhump Phase into Orbital Light Curve 

To proceed further, we need to make some more simpli- 
fying assumptions concerning the spatial distribution of 
the superhump light. We here assume that the disk rim is 
a circular cylindrical surface with radius R, and vertical 
width 77, and that the superhump light distribution at the 
disk rim is now represented by 1 = I{0), where / is spe- 
cific intensity of radiation and is the azimuthal angle of a 
surface element of the vertical rim in the corotating frame 
of the binary, respectively. Here, the azimuthal angle is 
measured in the direction of the binary orbital motion in 
our convention. That is, all of non-axisymmetric effects 
of the superhump light are assumed to be represented by 
an azimuthal dependence of the intensity of radiation at 
the circular disk rim. 

This assumption is admittedly very crude as the disk 
rim is not certainly circular, but is elongated during the 
superhump era. Here, we have restricted our attention 
only to the azimuthal distribution of the superhump light 
by neglecting all other effects. That is the reason why I 
have adopted the above-mentioned assumption. Because 
of this shortcoming, we must be careful not to over- 
interpret the obtained results. 

We represent the direction to the observer by {i,0o), 
where i is the orbital inclination and 0q is the azimuthal 
angle to the observer in the corotating frame. We here 
take into account the limb-darkening elfect of a surface 
clement, which is written by 



where /c is the intensity of radiation along the surface 
normal, •& is the angle between the surface normal and 
the direction to the observer, and P is the limb-darkening 
coefficient. The angle between the surface normal and the 
direction to the observer, is expressed by 

cos t9 = sin i sin 6*' , (2) 

where 0' is an azimuthal angle of a disk surface element in 
the observer's frame in which the origin is chosen to the 
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direction perpendicular to the observer. The azimuthal 
angles in these two frames are related to each other by 

0' = e-0o + ^. (3) 

Let us consider some certain supcrhump light distribu- 
tion from the disk rim at a given superhump phase, which 
is represented by 1^ = Ic{Q)- We may expand it into its 
Fourier components as 



(4) 



where m = 0, 1,2,3, .... are zero and positive integers and 
(/),„ is a constant phase for a given m. 

In what follows, we discuss each of Fourier compo- 
nents separately. The radiation emitted to the ob- 
server from length Rd9 of the cylindrical surface is 
RH cos -dim cos{m9 + <j>m)d9 for the m-th Fourier compo- 
nent. We now integrate the light contribution for the visi- 
ble surface of the cylinder, for which 9' varies from to it. 
The m-th component of radiation received by the observer 
will then be 



F„ 



l + f3 



X / {l + (3cosi})cos'dcos{m9 + (j>,n)d9', (5) 
Jo 

where D is the distance of the observer from the star. 
Substituting equations (2) and (3) into equation (5), we 
obtain 

RHI,JD^ 



F, 



1 + p 

X / (sinisin6'' + /3sin^isin'^6'') 
Jo 

mO' + m9o + 0m —j dO'. 

After some elementary calculations, wc can show that 
equation (6) can be written as 

RHIrrjD^ . 



(6) 



Fm{do) = 

x{Ajn+ PsmiBr,i)cos{ineo + 4im), 



(7) 



where Am and Bm are numerical constants. 

Equation (7) is the very equation which wc have sought, 
and it tells us that the orbital light variation due to the 
m-th component of superhump light is proportional to 
the intensity at the very surface element of the disk rim 
facing to the observer, i.e., the surface clement at the az- 
imuthal angle (?o- The coefficients A,„ and represents 
how much reduction of light occurs due to integration over 
the visible surface of the disk rim. As shown below, they 
decrease rapidly by going to the higher order m-mode be- 
cause of cancellation effects with positive and negative 
light contributions. 

It is found that the coefficient A,„ is given for an even 
integer m by 



Table 1. CocfBcicnts Am and Bm, and coefficient Cm for 
two cases with no limb darkening and with Umb darkening of 
13 = 0.56 and i = 90°. 



m 


Am 




Cm 


Cm 








(/3 = 0) 


(/3 = 0.56) 





2 


7r/2 


1 


0.923 


1 


7r/2 


4/3 


0.785 


0.743 


2 


2/3 


7r/4 


0.333 


0.355 


3 





4/15 





0.048 


4 


-2/15 





-0.067 


-0.043 


5 





-4/105 





-0.00068 


6 


2/35 





0.029 


0.018 



Am = (-1 



|l+m/2_ 



1' 



(8) 



while Am = for odd integers m, except for m = 1. The 
case for to = 1 is an exception and Ai = tt/2. Similarly, 
the coefficient Bm is given for odd integers by 

i3m = (-l)(-»+i)/2 l (9) 
mym'^ — 4j 

while Bm = for even integer to, except for to = and 
TO = 2, for which Bq = 7r/2 and B2 = 7r/4. Equation (7) 
may be written in a more convenient form as 

2RH 



Fm{do) 



= cos{m9o 
where the coefficient 

Am+ PsiniBm 



sini/„iC„i cos{m9o + 4>m) 



(10) 
(11) 

2(1 + /?) 

is the transformation coefficient from the brightness dis- 
tribution of the disk rim to the orbital light curve. It is 
normalized to one for a disk rim with uniform brightness 
and no limb darkening. Their numerical values arc listed 
in table 1 for two cases with no limb darkening and with 
limb darkening of /3 = 0.56 and an inclination oi i = 90° 
together with coefficients Am and Bm- 

3.3. Inversion: Reconstruction of Azimuthal Brightness 
Distribution of the Disk from Orbital Light Curve 

The inversion from the orbital light curve to the az- 
imuthal brightness distribution is straightforward in our 
case. Let us write an orbital light curve as = i^(<po) 
for a certain superhump phase obtained by the method 
described in section 2, where tp^ is the orbital phase in 
angular units. The orbital phase tp^ and the azimuthal 
angle of the observer in the corotating frame of the bi- 
nary, 9q, used in the previous subsection arc related each 
other by (/jq = 27r — 9q. Since F — F{lpq) is a periodic func- 
tion with 27r, we may expand the orbital light curve into 
Fourier series as 



F^^ FmCos.{mipo - 4>r, 



(13) 



m=0 
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Table 2. Observed power spectra of super hump light curves 
in the first week of the common super hump era for the 2001 
outburst of WZ Sge (data from Patterson et al 2002). 



n Frequency Power Relative amp. Identification 

(c d~^) On 

1 17.49 1030 32 ujQ-n 

2 17.64 ^ 60 TVT ujq 

3 34.97 ~ 85 9.2 2ljq - 2n 

4 35.12 ~ 35 5J 2tJo - O 

5 52?77 ^^50 7A 3uja-n 
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there is no ambiguity of tlieir identification. Two largest 
peaks atLUo-n= 17.49 c d"^ and 2uJo - 20 = 34.97 c d"^ 
are those of the superhump period and its first harmon- 
ics, corresponding to the superhump light curve shown by 
Patterson et al.(2002) at the lowest left frame of their fig- 
ure 6, which was obtained by synchronous summation at 
the superhump period. The most interesting is the signal 
at 3ujQ — ft = 52.77 c d~^ which appears as an isolated 
peak at the frequency range of the second harmonics and 
its side bands. As will be shown below, it is identified a 
traveling two-armed spiral wave. 

The light intensity in the first week of the common su- 
perhump era is expressed as 

m 

— Lo + ai cos [{ujQ - 0)t + tpi] + 02 cos(u;oi + V'2) 
+03 cos [2{ujo — + V'3] + 04 cos [(2ti>o — + i^i] 
+ar,cos[{3ujQ-n)t + 'ip5], (16) 

where Lq stands for the background light and and ipi 
are the amplitudes and constant phases for the respective 
modes. In equation (16) and following mathematical ex- 
pressions, we understand that the frequencies, such as luq 
and S7, and phases, such as ipsa and ipheat, are given in 
angular units because there may be no confusion, for in- 
stance, with cyclic frequencies. Equation (16) is rewritten 
for later convenience as 

Lit) 

— Lq + qi cos [{uiQ — fl)t + ipi] 

+a3 cos [2(^0 - 0)t -I- ■03] + 02 cos(woi + ^2) 
4-04 cos [ujot + (luq — fl)t + Ipi] 

+a5Cos[2ujQt+ {LJo-n)t + tp5]. (17) 

Let N denote the number of the orbital light curves in 
the beat cycle, which is given by the nearest integer of the 
ratio ujo/n. A particular time in the n-th light curve is 
then expressed as 

u;ot = 2'K{n-l) + tpo, (18) 

{uo-n)t^2TT{n-l) + ipsii. (19) 

where ipo and (psn are the orbital phase and the super- 
hump phase, respectively. They are related to each other 
by ipo = ipsa + (/^boat, where t/Sboat = is the beat phase. 
Substituting equations (18) and (19) into equation (17), 
we obtain 

L = Lo + ai cos((y9sH + V'l) + 0,3 cos(2(y9sH + ips) 
+a2 cos((^o + V'2) + 04 C0s((y9o + ^SH + V'4) 
-ha5COs(2(/3o + </'SH + V's)- (20) 

Let us now use the concept of helical tomography, in 
which each "orbital light curve" corresponds to a fictitious 
light curve of a snap shot of the accretion disk [see, e.g., 
ten snap shots of the accretion disk in figure 7 of Hirose 
and Osaki (1990)] viewed from different azimuthal direc- 
tions. We may regard equation (20) as that describing an 
"orbital light curve" at a given superhump phase. That 
is, in the above equation ipsn is regarded to be a fixed 



(15) 



= ^ F^cos{m0o + <f>m), (14) 

where F^^ is a Fourier amplitude with m-th component. 
Here is a constant phase of the m-th component and 
a minus sign is attached to 4'm in equation (13) because 
it has turned out that the constant phase 4>m defined in 
this subsection is identical to that defined in the preceding 
subsection. We note here that the zeroth component, Fq, 
includes not only contributions from the vertical disk rim, 
but also from the horizontal disk surface. 

We now recover individual mode amplitudes, Im, of the 
azimuthal brightness distribution of the disk by inverting 
equation (11) to obtain 

/ -( 
" \2RHsmiJCm' 

Once we obtain all Fourier amplitudes, we reconstruct the 
disk image by adding all Fourier components in equation 
(4). 

4. New Light on the Power Spectra of 
Superhump Light Curves for the 2001 
Outburst of WZ Sge 

Patterson et al. (2002) made a power spectral analysis 
of superhump light curves of the 2001 outburst of WZ Sge 
over the full beat cycle, and obtained many frequencies; all 
were linear combinations of coq and 0, where Uq and O are 
the orbital frequency and the putative apsidal frequency 
of an eccentric disk, respectively. As a matter of fact, we 
do not even need to construct new "orbital light curves" 
to demonstrate the power of the helical tomography if we 
already have information of power spectra of superhump 
light curves. 

We explain this by using Patterson et al's (2002) data 
for the first week of the common superhump era of the 
2001 outburst of WZ Sge. Table 2 summarizes results of 
power spectral analysis taken from their figure 6. Each 
column of table 2 indicates from left to right, the ordinal 
number, n, of individual modes with increasing frequency, 
cyclic frequency in units of cycle d~^, power of modes es- 
timated by eye from their figure 6, relative amplitudes of 
modes, o„, and its identification. As can be seen in their 
figure 6, five discrete peaks in power spectra are seen and 
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constant, while (/jq = 27r — 9q is a variable orbital phase, 
representing the azimuthal angle of the line of sight of an 
observer in the co-rotating frame of the binary. 
We may rewrite equation (20) as 



L = Lq + CIi COs((psH + 1pl) + a3 C0S(2(^SH + "03 ) 

-|-a2Cos(6'o-V'2) + a4Cos(6'o -ipsu^^Ji) 
4-05 cos(26'o - ipsa - V's)- 



(21) 



By using the method described in the preceding section, 
we now invert equation (21) to obtain the surface bright- 
ness distribution of the vertical disk rim as 



i"(6',VSH 



/2RHsmi 



Lo + bi C0S(</7SH + 0i) + ^3 cos(2v?sH + V'3) 
-1-62 cos(6' - 02) + 64 cos(6' - ipsn - 04) 

-^65 COs(26'-V3sH-05), 



(22) 



where 61 = ai, 63 = 03, 62 = ailC\, 64 = a^jCx, and 65 = 
o-^jCi in the case of no limb darkening. 

Equation (22) represents the azimuthal brightness dis- 
tribution of the superhump light at the disk rim for a given 
superhump phase. In order to reconstruct the brightness 
distribution of the disk for a certain superhump phase, 
we need information about the constant phases, 0„, of 
individual Fourier components. Unfortunately, this infor- 
mation is usually not available in published form. Since 
the phase information of individual modes has turned out 
to be very important from the standpoint of the helical to- 
mography, we urge observers to publish such information 
together with the amplitudes of individual modes. 

Since v^sh = (wo — ri)i, we may rewrite the above equa- 
tion as 



/(0,t)( 



/ 2RHsmi 



1)2 



= Lg + bi COS [(^0 — + -01] + ^3 cos [2{ujo ^ + 03] 
+b2 cos(6'-02) 
+b4 cos [9 — (cjo — — 04] 

+b5COs[29-{uja-n)t~ip5]. (23) 

Equation (23) describes the brightness distribution of the 
superhump light as a function of the azimuthal angle, 9, in 
the disk and time, t. As expected, it is a periodic function 
with the superhump period. 

The first line of the right hand side of equation (23) 
is then a certain constant at a given superhump phase, 
representing the light level of a synchronously summed 
superhump light curve (see, the lowest left panel of figure 
6 of Patterson et al. 2002). The second line of equation 
(23) represents a one-armed wave fixed in the corotating 
frame of the binary. It corresponds to the lowest right 
panel of the same figure; however a much clearer figure 
is found in the lowest right panel of figure 7 for the sec- 
ond week of the common superhump era (Patterson et al. 
2002). This was called "mean orbital light curve", and 
was interpreted as a hot spot light curve by Patterson et 
al. (2002). As discussed by Osaki and Meyer (2003), we 
have argued that the "mean orbital light curve" simply 



represents the non-axisymmetric distribution of the su- 
perhump light, and that it is not due to the enhanced hot 
spot. This point will be taken up once more in the next 
section. 

The third line also represents a one-armed wave, but 
this time it is a traveling wave which rotates progradely 
once per superhump period in the co-rotating frame, be- 
cause 

d9 , , 

-..o-f^. (24) 

If we introduce the azimuthal angle in the inertial frame 
by 8, it is related to that in the corotating frame by O = 
9 + LOQt. Thus, this wave mode rotates twice in the inertial 
frame per superhump period. 

From the same discussion as given above, we find that 
the last term of the right-hand side represents a two-armed 
traveling wave, as 

d9 _ ujQ-n 

That is, the wave pattern rotates one half per superhump 
period in the corotating frame. However, the same disk 
configuration returns back after one superhump period 
because the basic pattern is two-armed. 

The brightness distribution of the disk can also be ex- 
pressed in the inertial frame (i.e., in the observer's frame). 
Since the azimuthal angle in the inertial frame, 8, is re- 
lated to that in the corotating frame. 9, by = 8 — ujot, 
equation (23) is rewritten as 

, f2RHsmi 

= Lq + feicos[-(a;o - - 0i] + &3cos[-2(ti;o - - 0'3] 

+b2 C0S(8 - UJot - 02) 

-I-64 cos [8 — {2uJo — i})t~ 04] 

cos [28 - (3^0 - f^)i - 05] • (26) 

Equation (26) describes the brightness distribution of the 
disk in the inertial frame. 

The above results are summarized in table 3 where wave 
modes are specified by {m,£) in which the brightness dis- 
tribution in the disk is written as cos(m8 — tujot)- Here, 
a small frequency shift of f2 to iujQ is disregarded in speci- 
fying the wave modes, because n ^ujq. The mode identi- 
fication of five modes with (m,£) is straight forward from 
equation (26). The main peak of the power spectrum (i.e., 
the n = I mode) is the well-known superhump periodic- 
ity with frequency loq — fl and its wave is identified by 
the (O,l)-mode because there is no dependence on the az- 
imuthal angle, 8, as seen from equation (26). The slowly 
precessing eccentric disk mode with cos(8 — ilt)-wave is 
identified as the (l,0)-modc but this mode, itself, does 
not produce any brightness distribution except for secu- 
lar variation with the beat period (i.e., 27r/il). The tidal 
coupling of the eccentric (l,0)-mode with the compan- 
ion's gravitational field described by (m, TO)-waves pro- 
duces (m — 1, 7Ti)-modcs as discussed in Lubow (1991a) 
and in the next section. The superhump light variation 



No. ] 

Table 3. Mode identifications and inverted amplitudes of in- 
dividual modes. 



n Frequency Identification Mode Inverted amp. 

(c d-i) (m,e) bn 

1 17.49 LOo-n (0,1) 32 

2 17.64 luq (1^1) 9.8 

3 34.97 2LOo-2n (0,2) 9.2 

4 35.12 2u>o-n (1,2) 7.5 

5 52.77 Sluq - n (2,3) 22 



with the (0, l)-mode is produced by a tidal interaction be- 
tween the eccentric (l,0)-mode and the tidal (l,l)-mode. 
The various wave modes excited in the superhump disk 
are discussed in the next section. 

We find from table 3 that the wave mode with cos(20 — 
3uj()t) was very strongly excited in the first week of the 
common superhump era in the 2001 outburst of WZ Sge; 
its significance is discussed in the next section. 

The above result can easily be extended to a Fourier 
component with frequency hloq — fl. Since the time de- 
pendence of this mode is written as 

cos [{nujQ - n)t + %jjn] = cos [{n - l)ipo + fSH + V^n)] ,(27) 

it represents an (n — l)-armed traveling wave which rotates 
[l/(n — l)]-times per superhump period. However, the ba- 
sic disk configuration returns back after one superhump 
period because of its (n — l)-armed structure. Patterson 
et al. (2002) found in the second week of the common su- 
perhump era of the 2001 outburst of WZ Sge the primary 
Fourier components at nujo — f2 in the power spectrum 
with n = 1, 2, 3, 4, 5, 6, 7, 8, 9. These are exactly (n — 1)- 
armed waves discussed here. The appearance of higher 
Fourier components in the second week indicates the ap- 
pearance of a much shaper feature in the second week for 
the azimuthal brightness distribution in the disk. 

5. Lubow's Theory of Non-Linear Wave Coupling 
by the Eccentric Lindblad Resonance 

In the previous section we have demonstrated that by 
using the concept of helical tomography we can identify 
azimuthal wave numbers of various frequencies found in 
the power spectral analysis of superhump light curves by 
Patterson et al (2002). In particular, we identify the 
Fourier components at nujQ — D, as {n — l)-armed trav- 
eling waves. In this section, we discuss why these wave 
modes are excited in the accretion disk of a superhump 
binary. To do so, we use Lubow's theory (1991a, b) of non- 
linear wave coupling. Below, we basically follow Lubow's 
(1991a, b) discussions on this problem. 

In this section, we use mostly the inertial frame of ref- 
erence rather than the corotating frame. The tidal per- 
turbing potential in the accretion disk by the secondary 
star is expressed as 

(f>{r,Q) ^ (t)r,i{r) cos{m9) = 0,„(r) cos [m(e - wq^)] . (28) 
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Let us consider wave fields in the accretion disk which are 
specified by two integers (m,£) in which its azimuthal and 
time dependence is expressed as cos {m& — ituot) . The 
tidal perturbation potential produces forced wave fields 
with m ^ m and £ = m, as is known very well. 

Let us now impose some finite eccentricity in the ac- 
cretion disk. The disk eccentricity is described as a one- 
armed eigenmode whose eigenfrequency, fl, is known to 
be very small as compared with ljq (see, Hirose, Osaki 
1993). Following Lubow (1991a), we denote this mode as 
the (l,0)-mode. The non-linear wave coupling between the 
(l,0)-modc and the (m,m)-mode produces another wave 
mode with (m — 1,to). Lubow (1991a,b) has demonstrated 
that, in particular, a strong wave with (2,3) is launched 
at the 3:1 resonance region of the accretion disk by the 
eccentric Lindblad resonance. Further coupling between 
the (m — l,m)-mode and the (m,m)-mode reinforces the 
imposed eccentric mode with (1,0). This is Lubow's pic- 
ture of excitation of the eccentric mode by the eccentric 
Lindblad resonance. 

So far we have neglected the small, but finite, frequency, 
fl, of the eccentric mode to describe various (m,£)-modes. 
Let us now take into account slow precession of an ec- 
centric mode (1,0). This mode has a small, but finite, 
frequency of fl, describing a beat period of about 6 days 
in the case of WZ Sge. Since the (to — l,TO)-mode is pro- 
duced by coupling between the (l,0)-mode and the (to,to)- 
mode, its frequency is given by towq — ^- This is exactly 
the (to — l)-armed traveling wave discussed in the previous 
section. 

The (0, l)-mode with frequency luq — ^ corresponds to 
the main peak of the power spectrum of the common su- 
perhump era with the superhump period, only a period- 
icity to be observed in superhump binary systems with a 
low orbital inclination. Many peaks in the power spectra 
found by Patterson et al (2002) in the common super- 
hump era are those with towq — which are interpreted 
as the (to — l,m) mode in the previous section. The dom- 
inance of these modes in the common superhump era is a 
natural consequence of Lubow's (1991) theory. In partic- 
ular, the fact that the (2,3) mode at a frequency 3a;o — ^ 
is very strong [see, the right most panel of the middle 
row in fig. 6 of Patterson et al. (2002)] in the first week 
of the common superhump era, supports Lubow's theory 
because the excitation of the (2,3)-mode plays the most 
important role for the growth of the eccentric mode by 
the 3:1 resonance in his theory. The (2,3)-modc is two- 
armed spiral waves, rotating one and half times around 
the accretion disk over one superhump period in the in- 
ertial frame, as clearly seen in hydrodynamic simulations 
by Simpson and Wood (1998). 

Another interesting aspect is a strong excitation of a 
mode with frequency ujq, in particular, in the second week 
of the common superhump era of the 2001 outburst of 
WZ Sge. This is clearly seen in the lowest right panel 
of figure 7 by Patterson et al. (2002), and it was called 
the "orbital light curve" because it was obtained by syn- 
chronous summation at the orbital period. Patterson et 
al (2002) interpreted it as being a hot-spot hght curve. 



Helical Tomography for Superhump Light Curves 



8 



Y. Osaki 



and they used it as evidence for enhanced mass transfer. 
This mode is a one-armed wave pattern with (1,1) fixed 
in the corotating frame of the binary, as discussed in the 
previous section. 

In Lubow's theory for the excitation of eccentricity of 
a disk, it is a two-step process. In the first stage, the 
(l,0)-mode grows rapidly together with the (2,3)-mode 
launched by the 3:1 resonance. In the second stage, once 
the (l,0)-mode attains a large amplitude, various non- 
linear wave couplings become important, and a kind of 
steady state is established between various wave modes. 
In this picture, the (l,l)-mode is also reinforced by non- 
linear wave coupling between the original eccentric (1,0)- 
mode and the (O,l)-mode responsible for the superhump 
light variation. The (l,l)-mode observed with a large am- 
plitude in the second week is thus a result of superhump 
phenomenon itself. If so, the amplitude of the (l,l)-modc 
will vary in a long run together with the amplitude of 
the (O,l)-modc; in other words, the amplitude of "orbital 
wave" in Patterson et al. (2001) should vary together 
with the superhump light amplitude, a picture already 
presented by Osaki and Meyer (2003). 

6. Application to Other Systems 

An SU UMa-type dwarf nova, lY UMa, is another 
eclipsing binary system for which Patterson et al. (2000) 
found a rich spectrum of frequencies in the power spectral 
analysis during its superout burst. Its power spectrum is 
very similar to that of WZ Sge, because signals at mujo — ^ 
with m = 3,4,5 were found besides the well-known super- 
hump periodicity at — i^. They are identified as (m — 1)- 
armcd traveling waves, just in the same way as in the case 
of WZ Sge, as discussed in sections 4 and 5. 

There exist several eclipsing SU UMa stars. The most 
well known among them are OY Car, Z Cha, and HT Cas. 
Their superhump light curves could be examined by new 
light based on the concept of helical tomography. 

The AM CVn stars are thought to be cataclysmic vari- 
able stars having helium-rich and hydrogen-deficient ac- 
cretion disks (Warner 1995). In particular, its prototype 
star, AM CVn, itself exhibits very complex light varia- 
tions. As already pointed out by Solheim and Provencal 
(1998), the disk structure will be analyzed by examining 
a rich spectrum of light curves of AM CVn stars. By 
analyzing its light curve, Skillman et al. (1999) have ob- 
tained periodic signals as many as 20 discrete frequencies. 
According to Skillman et al. (1999), the basic clocks of 
this star arc an orbital period of 1028 s and a superhump 
period of 1051 s. However, many of the discovered fre- 
quencies can be understood as harmonics and sidebands of 
these fundamental frequencies, i.e., with hloq — mfl, where 
ujQ and f2 are the orbital frequency and the apsidal preces- 
sion frequency of an eccentric disk, respectively. On the 
other hand, Solheim ct al. (1998) have made a quite dif- 
ferent interpretation in that the fundamental frequency at 
951/U,Hz (with a period of 1051s) is related to the orbital 
period of the binary and the dominant light signal at 525 
s is due to the two-armed spiral structure in the disk fixed 



in the corotating frame of the binary. 
7. Summary and Discussions 

We have proposed a new method to analyze complex 
superhump light curves of SU UMa stars with high or- 
bital inclination, particularly for the 2001 outburst of WZ 
Sge. It is called the helical tomography because it can 
reconstruct a series of accretion disk images in various 
superhump phases. Based on the concept of helical to- 
mography, power spectral data of superhump light curves 
obtained for the 2001 outburst of WZ Sge by Patterson et 
al. (2002) have been interpreted, and the azimuthal wave 
number of various frequency modes identified. 

As already discussed in section 3, a very crude assump- 
tion had to be made for the superhump light distribution 
in order to perform an inversion of the orbital light curves 
into the azimuthal distributions of the superhump light. 
This is one of the biggest shortcomings in our theory and 
a much better method should be explored. Nevertheless, 
we can demonstrate that a two-armed traveling wave with 
cos(20 — 3ujot) was excited vigorously in the first week of 
the common superhump era, supporting Lubow's (1991) 
theory of the eccentric Lindblad resonance for the tidal 
instability. 
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